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Abstract 

We have developed a theoretical method for interpretation of photoionization experiments with 
the H3 molecule. In the present study we give a detailed description of the method, which combines 
multichannel quantum defect theory, the adiabatic hyperspherical approach, and the techniques of 
outgoing Siegert pseudostates. The present method accounts for vibrational and rotation excita- 
tions of the molecule, deals with all symmetry restrictions imposed by the geometry of the molecule, 
including vibrational, rotational, electronic and nuclear spin symmetries. The method was recently 
applied to treat dissociative recombination of the HJ^ ion. Since dissociative recombination has 
been a controversial problem, the present study also allows us to test the method on the process 
of photoionization, which is understood better than dissociative recombination. Good agreement 
with two photoionization experiments is obtained. 
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I. INTRODUCTION 



The simplest polyatomic molecules, H3 and H3 , have been intensively studied during the 
last decades. Interest in these molecules is motivated by the fact that the ion plays an 
important role in the chain of chemical reactions in interstellar space, acting as protonator 
in chemical reactions with almost all atoms. In particular, dissociation recombination (DR) 
of with an electron leads via several intermediate steps to the production of water in 
interstellar space. Many successful models in interstellar chemistry are based on Hjj' DR. 
In addition, H;^ attracts theorists as a benchmark for high accuracy calculations with small 
molecules: Theoretical ab initio methods can be tested against existing experimental H3" 
spectroscopy data. The interest in the neutral H3 molecule is closely related to the problem of 
Ho" DR. But H3 also presents great interest from another point of view. Experimental studies 
U 0, 0, y, 0| of the metastable H3 molecule and several later theoretical studies ^ Is, 0, [t] 
revealed non-Born- Oppenheimer effects of coupling between its electronic, vibrational, and 
rotational degrees of freedom. As was shown recently 0,0], these non-Born-Oppenheimer 
effects play an important role in H3" DR as well. The present study is devoted to a theoretical 
treatment of II3 photoionization. 

There are three main reasons for this study. The first one is to explore a new theo- 
retical method for the treatment of polyatomic photoionization. Our method is based on 



multi-channel quantum defect theory (MQDT) 
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12l IS^ . the adiabatic hyperspherical 



approach to vibrational dynamics of three nuclei, the formalism of outgoing wave Siegert 
states and inclusion of a non-Born-Oppenheimer coupling — Jahn- Teller effect. 



The second reason is related to a recent study of Ht DR 



where the reported method 



was very successful in treatment of H;^ DR, giving good agreement between theoretical 



d a gree ment 



However, since that 



calculations and experimental results from storage rings 
method is new, it is desirable to test it in greater detail. An application to the interpretation 
of H3 photoionization experiments k , ^ is such a test. These photoionization experiments 

nnriQ 

were successfully interpreted in previous theoretical work |^ 13, M, 12I| , where another method 
based on MQDT was applied. Thus, the present treatment can also be tested against the 
previous theoretical studies. 

Our treatment of photoionization is similar to the one developed by Stephens and Greene 



and employed in Refs. 



for interpretation of two photoionization experiments 



2 



full rovibronic frame transformation 



by Bordas et al. [4] and by Mistn'k et al. p. Both experiments were interpreted using a 



The present treatment has several differences 



from the one proposed by Stephens and Greene. The first difference is the use of the 



adiabatic hyperspherical approximation 



2G 
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22| | for the representation of vibrational 



wave functions. Stephens and Greene used the exact three-dimensional vibrational wave 
functions. The second difference is the correction of the incompatibility between the form 



for the reaction matrices used in Refs. 



23l | and the quantum defect parameters 



of Jahn- Teller coupling used in the mentioned studies. In fact, the values of Jahn- Teller 
quantum defect parameters used in Refs. compatible with an alternative 

form of the reaction matrix, which was adopted in Refs. [2J, |25|, |26|. In the present work 
we use the same form of /T-matrix as in Refs. 0, 0, Q, |2^ and quantum defect parameters 
from Ref. P]. Thus, Jahn- Teller parameters 5 and A from should be multiplied by — tt 
to be used in the present study. The third difference is in the symmetrization of the total 



rovibrational wave functions of the Hg ion. In Refs. 



the symmetrization is made 



according to the procedure proposed by Spirko and Jensen j22|: Rotational and vibrational 
parts of the total wave function are symmetrized separately and in two-step procedure. In 
the present treatment we symmetrize the total wave function only once at the very final 
step. This greatly simplifies the construction of wave functions of a required symmetry. 
The fourth difference is in calculation of dipole transition moments. Calculating the dipole 
moment into a final state, Stephens and Greene accounted only for the diagonal component 
of the final state wave function. Our treatment accounts for all non-diagonal wave function 
components contributing to the dipole transition element. 

The article is organized as follows. Section |n] describes construction of the total wave 
function of Hjj' and compares our method of the construction with the method proposed in 
Ref. {2^. In Sec. IIIH we build up the scattering matrix that represents the collision between 
an electron and an ion. Section IIVI presents a derivation of dipole transition moments and 
oscillator strengths for H3. We discuss results of our calculation and compare those results 
with experimental data in Sec. |3 Section IVII states our conclusions. 

Atomic units are used in the article unless otherwise stated. 
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II. SYMMETRY OF THE TOTAL WAVE FUNCTION OF H+ 



In this study we consider only p-wave scattering (or half-scattering) of the electron from 
the molecule. As demonstrated in Refs. III higher electronic partial waves make much 
smaller contribution than the p-wave to the photoionization spectrum. Similar to our study 
of U.^ dissociative recombination, we chose the molecular axis Z along the main symmetr 
axis of the molecule. Directions of two other axes, X and Y, are shown in Fig. 3 of Ref. 



A. Total wave function 

The total wave function <|)"-*^™ of the ion can be represented as a sum of terms, each of 



which is product of three factors j9||: 

^n.sym ^ $^^7^^^^^^^ (c,, 7) $v( Q) • (l) 

In the above equation, a,f3, and 7 are three Euler angles defining the orientation of the 
molecular fixed axis with respect to the space fixed coordinates system. Below, we describe 
briefly the construction of all three factors in the product of Eq. (Q. A more detailed 
description is given in Ref. j9|. 

The rotational part n(a,P,-,) ot the total wave function in Eq. OJ^ the symmetric 
top wave function for Hjj', which is proportional to the Wigner function 28]. The quantum 
numbers A^"*", K~^, and refer to the total angular momentum A^"^ and its projections on 
the molecular Z-axis, K^, and the laboratory z-axis, m+. The transformation properties of 
the symmetric top wave function under the D^h group, are given in Table II of Ref. j^. 

The vibrational symmetry of B.^ and II3 is described by the group C^y. C^^ is a subgroup 
of D-^h- Dsh = o-h® Csj;, where ah is the operation of reflection with respect to the plane 
of three nuclei. For our discussion of the vibrational symmetry of H;^, it is convenient to 
use normal coordinates Qi, Qx-, and Qy (for definitions, see, for example, Ref. 0). Qi 
describes the symmetric stretch mode. The motion along this coordinate is characterized 
by the (approximate) quantum number Vi and by the corresponding frequency uJi. Normal 
coordinates and Qy correspond to two vibrational modes having the same frequency of 
oscillations uj2- Vibrations along and Qy are characterized by the approximate numbers 
Vx and Vy, correspondingly. The total vibrational energy can be approximated E = uji{vi + 
l/2)+uj2{vx+Vy+l). (The vibrational quantum numbers vi, Vx and Vy and the corresponding 



energy E are not exact as long as the ionic molecular potential is not exactly harmonic.) Due 
to the degeneracy of the Qx and Qy modes, the two-dimensional vibrational motion along 
the Qx and Qy coordinates can be equivalently represented in polar vibrational coordinates 
p and (f). Then, instead of quantum numbers and Vy, it is convenient to define V2 = v^ + Vy 
and I2, where I2 is associated with the motion along coordinate: vibrational angular motion 
around the symmetry axis. Thus, the vibrational energy is determined only by the quantum 
numbers Vi and V2'- E = uji{vi + 1/2) + uj2{v2 + 1). The number V2 shows how many 
vibrational quanta are in the asymmetric mode. The number I2 determines how many of the 
asymmetric quanta f 2 contribute to the vibrational angular momentum, —V2 < ^2 < "^2 • In 
reality, due to the anharmonicity of potentials, the vibrational energy of states with same vi 
and V2 but different I2 are slightly different. However, pairs of states with ±/2, where I2 7^ 3/c 
(here and below, k is any integer number), are strictly degenerate. This is a consequence 
of the fact that the D^h symmetry group has doubly degenerate representations. Thus, 
vibrational wave functions, 

$v(Q) = |t^i,4'), (2) 
of the ion are specified by the triad of quantum numbers f 1, v^2- The quantum number I2 can 
have values — f 2, — f 2 + 2 . . . f 2 — 2, f 2 and it controls the symmetry of the vibrational wave 
functions. States with I2 = 3k, with an integer k 0, can be of Ai or A2 symmetry. In order 
to distinguish the two symmetries using the number I2, we will label states Ai with positive 
I2, h = 3k, k > and A2 states with negative l2- A pair of states with I2 7^ 3k having both 
signs of I2 constitutes the degenerate pair of functions, that transform according to the E 
representation. In contrast to the numbers fi,f2^, the classifications with symmetry labels 
Ai, A2, or E are exact. 

In the present treatment the relative phase of degenerate states \vi,V2) with I2 7^ 3k is 
slightly different from the one used in Ref. Q]- Transformations of V2) with the present 
choice of phases are summarized in Table III of Ref. jol. 

The third factor of the total ionic wave function is the nuclear spin wave function. The 
nuclear spin states are classified according to the total spin I = 3/2 or J = 1/2. These 

n 

states are constructed as described in Ref. (2|. The result is $o(-^ — 3/2), = 1/2), 

and = 1/2) wave functions, transformed according to the Ai (B E representation 

of the respective symmetry group of three identical particle permutations. The state 
$q(M = 3/2) transforms according to the Ai representation; the states =1/2) and 
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= 1/2) transform according to the E representation. 
The transformation of the total wave function ^^'^^"^ is determined by the quantum 
numbers K~^, I2, and gj . The final step in the construction of the total wave function is an 
appropriate symmetrization of ^^■'^y"^. Since the total wave function should be antisymmetric 
with respect to (12) we determine $j: 




In the above equation S2 = 1 for all vibrational states excluding A2, for which S2 = — 1; 
h = —^'2 if ^2 = 3A; ± 1, and I2 = I2 if h = 3/c. The condition for antisymmetry with respect 
to (12) is specified explicitly. If the symmetrization is trivial, i.e. both terms in Eq. (jH)) are 
identical, then the wave function is 

$t = C"^(^+,^2,^?/). (4) 

It is only possible if gi = 0, = and I2 = 1'2- 

The fermionic nature of nuclei also requires that the wave function <l>t should be anti- 
symmetric with respect to operations of (13) and (23). It is only possible if $j transforms 
according to the or A/^ representations of the D^h{M) group. This condition can be 
written as G = 3k, where G = + I2 + gi- The determination of the total symme- 
try has one exception from the above rule. Namely, when the symmetrization is trivial: 
h = —l'2i = 0. For this case, the rovibrational part of the product (jH) has Ai or A2 
symmetry, thus, gi can only be 0. Finally, the overall parity of the total state, which is 
determined as transformational under the operation of total inversion E*, is determined by 
the number K~^: The parity is even if is even and the parity is odd if K'^ is odd. 

In this study we are primarily interested in the ortho-modification of the H3 molecule of 
A'2 symmetry. Thus, Eq. is reduced to 

$f = (iZN+K+m+ia, f3,l)\Vl,V2^) - (-1)^'^ S2\'RN+-K+m+{<^, P,'l)\Vl,V2)) (5) 

or, when = and I2 = I2, to 

^t = nN+K+m+ioi,p,'j)\vi,V2^). (6) 

Only states with G = + I2 = 3k and with even are allowed. Again there is an 
exception, when the symmetrization is trivial: When K'^ = and is even (rotational 
symmetry is A'^), I2 = 3k can only be negative {A2 vibrational symmetry); when = 
and A^"^ is odd, I2 = 3k must be positive or zero. 
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B. Ht vibrational dynamics in an adiabatic hyperspherical approach 



As in our study of H3 dissociative recombination j8|, |9|, |23|, we employ the adiabatic 
hyperspherical approach to describe the vibrational dynamics of B.f and H3 in three di- 
mensions. In this approach, three hyperspherical coordinates, the hyper-radius R and two 
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22]. In our 



hyperangles 9, y9, represent three vibrational degrees of freedom , ^ 
calculations we use accurate potential surfaces of from Refs. 

The hyperspherical coordinates used are symmetry-adapted coordinates: Each operation 
of the Cat, group-a permutation of instantaneous positions of the three nuclei-is described 
with an appropriate change in the hyperangle ip only. With this respect, the hyperspherical 
coordinates are ..U. Jt.e „o.n,a, cooJina.. of Hj 033. wh.e a„ operation. 
C^y involve the polar angle (p uniquely. For example, the effect of the (123) operation is a 
cyclic permutation of the three internuclear distances. In the hyperspherical coordinates, it is 
realized by adding the angle 27r/3 to as determined in Eqs. (23) of Ref. |9|. The operation 
(12) exchanges the internuclear distances r2 and r^. Equations (23) of Ref. show that 
this operation corresponds to a mirror reflection about the axis ip(i2) = 7r/2, 7r/2]. Only 
the angle (p is changed, into (p' = {12)ip = tc — ip. This operation is exhibited in Fig. ^ The 
operation (12) exchanges nuclei 1 and 2, transforming the triangle Ta into Tb. The figure 
also shows all three symmetry axes, corresponding to three binary permutations (12), (23), 
and (13). These symmetry properties of the hyperspherical vibrational coordinates simplify 
our treatment appreciably. 

In the adiabatic hyperspherical method we first solve the vibrational Schrodinger equa- 
tion at a fixed hyper-radius R obtaining a set of energies Ui and corresponding eigen- 
functions ^i{6,ip). Changing R, we obtain a set of adiabatic potential curves Ui{R) and 
adiabatic hyperspherical eigenstates ^i{6,(p; R). As mentioned above, each element from 
the symmetry group is represented by a corresponding transformation involving only 
the hyperangle (p. The hyper-radius is not involved in the C^^ operations. Thus, the vi- 
brational hyperspherical states ^i{9,ip;R) and curves Ui{R) can be classified according to 
irreducible representations of the group C^y. Namely, each state ^i{9, ip; R) and correspond- 
ing curve Ui{R) can be labeled by either the Ai, A2, or E irreducible representation. The 
representation E is two-dimensional. Thus, two degenerate i?-components $j(0, </?) will be 
labeled by Ea and Eh- Their linear combinations are also good vibrational eigenstates. 



7 



Several low ionic hyperspherical curves U^{R) are shown in Fig. |21 We specify the 
pair v^2 quantum numbers for first several states. As in the familiar Born-Oppenheimer 
approximation for diatomic molecules, curves of the same symmetry do not cross, whereas 
curves of different symmetries may cross. One can see from Fig. |21that low-lying potential 
curves with the same V2 are almost degenerate. This is because the anharmonicity is quite 
small for low states. However, at f 2 = 4 the potential curves with quantum numbers f 2^ = 4^ 
and = 4*^'^ are already significantly separated. 

Similar to Ref . l9| , from real- valued i?-states obtained after a diagonalization at fixed R, 

n 

in the space of the two hyperangles 9 and v?, we construct "helicity" E states: 

\v[-) = l^{Ea + ^E,)■ 

\vh = \v2'') = ^{Ea-iE,). (7) 

The sign of I2 in the above equation is chosen in such a way that (123) transforms the state 

Iff ^) as 

(123)14'^) = e^|4'^). (8) 

For example, if I/2I = 2, then I2 = —2 and I2 = 2 (see discussion in Sec. Ill C|) . Finally, 
we multiply all real-valued vibrational functions A2 by i in order to obtain a real reaction 
matrix K. 

Once the adiabatic hyperspherical potential curves U^{R) are determined, we calculate 
vibrational energies, by solving the adiabatic hyper-radial equation, 

' 1 

where {v,i} = {^1,^2^}. In solution of 
obey outgoing wave boundary conditions 
normalized as in Ref. j^. 

Inclusion of Siegert states into the treatment allows us to represent the dissociation of 



+ U^{R)\^UR)=E,,,^+,{R), (9) 

(EI), we seek solutions (Siegert states) that 
isl at a finite hyper-radius Rq and which are 



the neutral H3 formed during the collision between B.^ and the incident e . 



C. Comparison with an alternative symmetrization procedure 



An alternative procedure for the symmetrization of the total wave function was proposed 
by Spirko et al. |27|. Spirko et al. describe how rovibrational wave functions of different 
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D^h representations, F = A^, A2, E'^, and E^, are obtained from the products of rotational 
and vibrational parts. We use the symbol C to specify the parity of a state, where ( — or 
=". (In Ref. j^], the symbol f was used for this purpose.) In the approach of Ref. [27 1, 
the rotational states with F = A^, -E^, or E\^ are obtained by combining the symmetric 
top states lA^"*", i^"*", m"*") (Eqs. (59)-(61) of Ref. j^])- Then products of rotational and 
vibrational states are constructed. At the next step, these products are symmetrized again 
to give rovibrational states of good symmetry. This final step is quite laborious since one has 



to consider al 
(77) of Ref. 



would have to construct nuc 
(Eqs. (62)- (77) of Ref. 



possible combinations of different rotational and vibrational states (Eqs. (62)- 
27l|). To properly include all states of different nuclear spin symmetries, one 
ear spin states and then use a similar symmetrization procedure 
) one more time. In implementing the procedure of Spirko et 



221, 

al. for construction of the e~+Hjj" scattering matrix, we have found that it is difficult 
to obtain all states of the right symmetry. When carried out incorrectly, our scattering 
matrix displayed non-zero matrix elements between states of different symmetries, which 
of cause signals an error. For this reason, we think that our symmetrization procedure is 
advantageous for the present study, where we typically include and symmetrize hundreds 
of states. Our procedure involves only two simple steps and can be easily automated on a 
computer. First, the products of rotational, vibrational, and nuclear spin non-symmetrized 
states are constructed. For each product, the number G [or the total symmetry for the 
states that are trivially symmetrized, such as Eq. (jHl)] is determined. It describes behavior 
of the product under the symmetry operations (12) and (123). Then the symmetrization 
procedure (if it is not trivial) is accomplished by a single equation [see Eq. (jB []. 

To facilitate the comparison with the procedure, proposed by Spirko et al. 27|, we intro- 
duce a "helicity" pair of degenerate E states, which transform in a uniform way, irrespective 
of their nature: rotational, vibrational, or nuclear spin. We determine two degenerate states 
\Ej^) and \E^) by their symmetry properties 

{123)\E±) = e^'^\E^), 

(12)|E±) = |E^), (10) 

irrespective of coordinates: rotational, vibrational, electronic or nuclear spin. Using Ta- 
ble II of Ref. j^, it can be easily verified that rotational \E^'^) states are obtained from 
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|A^+, i^"*", m"*") according to 

\El^) = \N+, K+, m+), if = 3k + 1, 
1^5) = (-l)^+|Ar+,-ir+,m+), ifir+ = 3^ + l, (11) 
1^:;;^) = |A^+, m+), if = 3k + 2, 

1^:'^) = (-l)^^|Ar+,ir+,m+), ifK+ = 3~k + 2. 

In the above equation, A; is a non-negative integer. Vibrational {E"^) states are obtained 
using Eq. ©as 

\El) = \v^,vf'), if |/2| = 3fc + l, 

\El) = \vi,vf'), if I/2I =3fc + 2. (12) 



In addition to \E±) states we introduce another pair of states \Ea) and \Eb) as 

1 
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E±) = ^{\Ea)±t\E,)). (13) 



or 



\Ea) = ^i\E+) + \E^)), 

\E,) = - 1^-))- (14) 

The states \Ea) and \Ei,) are real. Using Eqs. pO|) and p4j). we obtain that 

il2)\E,) = \Ea), 
{12)\E,) = -\E,). (15) 

Using Eqs. (|TT1|) . (fT!?jl . and (fT^ and the fact that e^*^ = — 1±2^, we have for the operation 
(123) 

{m\E^) = -^-\E^)-y^\E,), 

{123)\E,) = y^\E^)-^-\E,). (16) 

At this stage, we can compare the present convention for \Ea) and \Eb) states with 
;he one from Ref. Comparing formulas (fT3j) and ()16p with Eqs. (51)-(54) of Ref. 

, we conclude that the two conventions for \Ea) and \Eh) vibrational states coincide. 
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A comparison for rotational \Ea) and states should account for a different choice of 
coordinate axes made in J27| and in the present study. In order to compare with the present 
work, axes x and y in j27l | must be exchanged. This will affect Eqs. (60)-(61) in 0|: 
The factor (—1)"' in Eqs. (60)-(61) at the second term of the equations must be omitted. 
After this modification, Eqs. (60) and (61) in describe exactly the same states as the 
rotational states determined in Eqs. (fTTj) and (fT^ of the present study. 

Now we can derive formulas (62)- (77) of 27 [for the product of rotational and vibrational 
states. Consider, for example, Eq. (75) of [27[ describing an overall state composed 
from \Ea) and \Eb) rotational and vibrational states. Using Eq. (|T^ and the fact that 
\E+) = \E_) ® \E_) and \E_) 



iV2 



we derive 
1 



E-)) = -^{\E1W_!<)-\EI)\e:;<)) 



i2V2 L 



i\E:) - ^\EJ:)) (\e:^^) - ^\El'^)) - he:) + zi^n) (ie:^"^) + w^^) 



E:Wt'^) + \e:)\e:''^)) . (17) 



iV2 

Comparing with Eq. (75) of [27[ we see that the \Eh) function of Eq. (75) differs from 
our state by an overall sign. The sign is important since it affects the result of the (123) 
operation. (The corresponding \Ea) state (Eq. (70) 27J) has the correct sign.) Note that 
Eq. (75) of 27|) also has a typographical error: Instead of symbol A2 in the first term, there 
should be \Ea). 

The rest of the formulas, (62)-(77) in 27|, can be derived in a similar way. 



III. THE SCATTERING MATRIX FOR EN ELECTRON COLLIDING WITH H+ 

A. Short-range scattering matrix of H^^ + e in presence of Siegert vibrational 
pseudostates 

Once vibrational Siegert pseudostates are calculated, the scattering matrix, describing the 
collision of the electron with the vibrating Hjj' ion, can be constructed. A vibrational frame 
transformation can be used to calculate the amplitude Si^f for the scattering from one 
vibrational state i' to another i. Here, the indices i and i' enumerate vibrational states and 
states of different projections A of the electron angular momentum. The vibrational part of 
the indices is represented by the triad {vi,V2}- The pair V2 represents one hyperspherical 
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curve U^{R); the index vi enumerates the Siegert pseudostates lying within that curve. 
Therefore the amphtude Si^i> for the process 

e-(A') + H3+K,t;?) ^ e-(A) + H+(t;i,4^) (18) 

is calculated in two steps, similar to the two-step calculation of vibrational energies. First, 
we determine i?-dependent amplitude 



(i?) = ((<l>,,,JSA;A'(Q) $.^,/0), (19) 

where the double brackets means an integration over hyperangles at constant hyper-radius 
R, and Q represents three internuclear coordinates. The scattering matrix Sa-a'{Q) includes 

ihe Jahn- Teller interaction and is calculated from the reaction matrix K as described in Ref. 

^ (see Eqs. (18)-(20) of Ref. Q). 

The scattering matrix S{Q) in Eq. (|19|) has indices A and A' and represents an amplitude 
for the process: 

e-{l = 1, A') + H+(Q) ^ e-(/ = 1, A) + H3+(Q). (20) 

Therefore, S'a^a'(2) represents the scattering amplitude when the electron scatters from 
one channel A to another A', while the nuclei do not have time to move. Equation (j2(jp 
describes the short-range B.^ + e~ collision in the clumped-nucleus approximation, where 
nuclear degrees of freedom are not yet coupled to the electronic degrees of freedom. 
The equation for the second step reads similarly as: 



where brackets ()^ means t 
implied surface term 



le integration in the sense of Siegert pseudostates, i.e. with an 
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16|: 



{i^vi,v2,hiR)\Sv2,i2Ay2A'^'(^)\^^'iy2A(^))s ^ 

/ '^V-i,V2,l2^-^) ' ^V2,l2,J^;vL,lL,\ 

'{R)-^,'^y^,i'^{R)dR+ (22) 

^0 

. '^Vl,V2,h{Rf)Sv2,l2A\v'2,l2,^'i^f)'^v[,V2^l2(^f) 

kvi,V2,l2 + kv[,v!2A 

When the integral in the above equation is evaluated, the usual complex conjugation of the 
bra wave function iIj^^^^^^i^^R) is omitted. The quantity ky^^v2,h ^ complex wave number 
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obtained from the complex energy E^^^^^^i^ of the corresponding Siegert state |fi,f2^) |9l Il6l|: 

E{v^,V2,l2} = ^W,i;2,/2}/(2/^) + ^^2,l2 (23) 

where Dy^^i^ is the dissociation hmit of the corresponding adiabatic hyperspherical curve 
Uv^^i^{R). In the present approach D^^^i^ is approximated by a value of U^^^i^i^R) at large 
hyper-radius i?/, D^^^i^ = U^^i^{Rf). 

Due to the presence of Siegert states with complex eigenenergies, this electron-ion scat- 
tering matrix is not unitary. The non-unitarity accounts for the fact that the electron can 
become stuck in the ion, leading to the dissociation of the system into neutral products. 



B. Rotational frame transformation and the final short-range scattering matrix 

The electron-ion scattering matrix Si^i', constructed above, does not account for the 
possibility of rotational excitation of the ion. If the H.^ ion is initially in one rotational state 
(N^' , K^'), a collision with the electron can scatter the rotational state into {N'^,K'^). 
Thus, an element iSj^j/ of the total scattering matrix S describes a transition from one 
rovibrational state i' = {v[,V2^}{N^' , K~^') to another i = {vi,V2}{N'^ , K~^). In indices i 
and i', we do not specify quantum numbers that are conserved during the collision. These 
quantum numbers are the total energy E, the total nuclear spin I of Hjj', the total angular 
momentum of the system and its projection m on the laboratory z-axis. (see Eq. fj26|) 
below) and, finally the total symmetry of the system: A2 or 

The change in the rotational excitation (N^' , K^') — » (N\ K~^) is taken into account 
using the rotational frame transformation approximation [l^ : Such a transition occurs 
mainly when the electron approaches close to the ion. Since a basis of rotational functions 
exists for which the short-range rotational Hamiltonian is diagonal, the transition ampli- 
tude for (A^+ ,K+) ^ {N+,K+) can be described by considering the coefficients that link 
the long-range quantum numbers {N~^,K^) with the short-range quanta. The short-range 
rotational states are specified by the projection A of the electronic angular momentum / on 
the ion-fixed Z-axis and by the projection K of the total angular momentum = A^"*" + / of 
the neutral molecule on the same Z-axis. A^ and / are conserved quantum numbers in both 
rotational bases. (This is one approximation of our treatment, because in reality /-changing 
collisions can occur with a small amplitude.) Below, we present a detailed description of the 
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rovibrational frame transformation for the e~+H3 system and specify all quantum numbers 
in both regions of interaction between e~ and Hjj'. 

At large electron-ion distances, the system is described by the electronic angular momen- 
tum I and its projection A on the laboratory z-axis, by the total ionic angular momentum A^, 
its projection m"*" on the laboratory 2;-axis and its projection on the molecular symmetry 
axis Z. Correspondingly, we represent the wave function of the e~+Hj]' system by a product 
of ^^■■'^y"^ and Yix{9,(p) (at this stage we do not specify electronic radial part of the total 
wave function): 

nN+m+K+{af3^)Yix{e,^)V{vi,v'i)<l>j. (24) 

The angles 6, (p are spherical angular coordinates of the electron in the laboratory coordinate 
system (LS). 

At short distances, the most appropriate molecular states, i.e., the states that almost 
diagonalize the Hamiltonian, are specified by the projection A of / on the molecular Z 
axis; by three internuclear coordinates Q since these remain approximately frozen during 
a single collision; by the total angular momentum of the system A^, including the electron 
momentum; and by its projections on the molecular Z axis, K, and on the laboratory z-axis, 
m. Thus, the total wave function at short distances is 

\N,K,m,A;v,,vi') = {-iy-''nN+m+K+{c^f3l)Yij,{e' ,^')\Q)<l>i. (25) 

The angles 6', io' determine the position of the electron relative to the molecular coordinate 

n 

system (MS). The transformation between the two wave functions is 9| 

I A^+, ir+; A^, m; v,,v',') = J2 ^f- wl^' ^' ^' ^' ^^^^2), (26) 

A 

which can be considered as if two angular momenta A^ and / with projections K and —A were 
added to give the momentum A^"*" with the projection = K — A. The quantum numbers 
N,l,m are not changed by the rotational, Eq. ()26|) . nor vibrational frame transformations, 
Eqs. (Uniand (jH))). Therefore, they are good quantum numbers at short and long distances, 
within the approximation of this study. 

Note that all three projections A, K, can be negative or positive. Equation (jSHI) differs, 
for example, from the one in Ref. jSJ] where all rotational functions are symmetrized with 
respect to different signs of projections. We keep both negative and positive projections 
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explicitly in order to symmetrize products of electronic, rotational, vibrational, and nuclear 
spin components of the total wave functions at the very final step. As was mentioned above, 
this simplifies the symmetrization procedure. 

The total short-range scattering matrix can now be constructed using the frame trans- 
formation techniques: When the electron is far from the ion, the interaction Hamiltonian is 
diagonal in the basis of the long-range wave functions; at short distances, short-range wave 
functions almost diagonalize the Hamiltonian. The short-range Hamiltonian is not exactly 
diagonal in the basis of states of Eq. ()25|). It has off-diagonal elements in A, owing to the 
Jahn- Teller coupling. The following selection rules can be formulated: (i) The Hamiltonian 
can only couple vibrational states of the same vibrational symmetry and the same value 
of A, or (a) A p-wave electron can couple the rovibrational channels according to the rule 
(A = l,/2 = —1) ^ (A' = — 1,^2 = !)• These selection rules insure that the total symmetry 
r of the system is conserved during the collision. 

The actual form of the coupling matrix Sa,a'{Q) is given in Refs. 0,y,Q- 

The final 

scattering matrix is represented as 



^{N,K,m,l,I,r) 
N+ ,K+ ,vi ■,N+' ,K+' ,v[ 2 



A,A 



N+ ,K+ 
L-A':N,K 



V{v[,v'^■Q)S^,^iQ)V{v^,v^i■Q)dQ 



'^l,-A;N,K- 



(27) 



The integral in the above equation is evaluated according to Eqs. and (j2H). The 
scattering matrix of Eq. (|27p is diagonal over quantum numbers A^, K, I, m, I and F. In the 
above equation F specifies the total molecular symmetry F = A2 or A2 of the considered 
state; I specifies the total spin 1/2 or 3/2. Therefore, photoionization oscillator strengths 
can be calculated separately for all possible values of these quantum numbers. 

In practice, we calculate S using non-symmetrized states of the type (0) and ()24|1 . The 
symmetrization procedure of Eq. ^ is then performed directly on the scattering matrix. 



Let the total dimension of the matrix S 



N+ ,K+,vi ; Af +' ,K+' ,v[ ,v'^^ 



be Nfot X Nfot. For the full 



specification of the scattering process, the photoionization threshold energies of rovibra- 
tional states i = {vi,V2}{N~^ .JC^); i = 1,2, ■ ■ ■ Ntot are needed. We use accurate energies 

35|, |36[. For some excited rovibrational levels, however, where no 



available in the literature 



data exist, the energies were calculated using the adiabatic hyperspherical and rigid- rotor 
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approximations. Another group of even higher energies are found to be complex, as expected 
for our Siegert pseudostate representation. 

Equation (|77|) gives the Ntot x Ntot scattering matrix describing e^+Hg colhsions. This 
matrix will be used to calculate the photoionization oscillator strengths. 



IV. OSCILLATOR STRENGTHS FOR THE INTERPRETATION OF H3 PHO- 
TOIONIZATION EXPERIMENTS 

We apply our method to describe two photoionization experiments involving the H3 



pjy our 

]^ 0] . aotn experiments nave pi 



molecule |5| . Both experiments 
multi-channel quantum defect theory 



laye _greviously been successfully interpreted using 
The present approach differs from previous 



theoretical studies. Two main differences are (1) the adiabatic hyperspherical method em- 
ployed for vibrational degrees of freedom and (2) an inclusion of the previously missed factor 
(— tt) in the Jahn- Teller parameters 5 and A. Below we give a detailed description of how 
the dipole transition moments and oscillator strengths are calculated. 

A. Vibrational wave function of the initial state 

In the experiment by Bordas et al. j^], the spectrum of the photoionization process, 

H3 3s A'2 (iV' = l,K' = 0, {00°}) + cji ^ H+ + e", (28) 

was measured. Here, uji is the photon energy. In the second experiment by Mistrik et al. 
P], an similar process was investigated. Mistrik et al. investigated photoionization starting 
from a different initial vibrational state, 

H3 3s A'2 (A^' = l,K' = 0, {10°}) + ^ H+ + e", (29) 

where UJ2 is the photon energy. In both processes indicated, the initial symmetry refers 
to the total molecular symmetry. The total molecular spin is / = 3/2 in both experiments. 
The total symmetry can be viewed as a direct product of the symmetry of the Yi^ ion 
with A^^ = 1, = 0, / = 3/2 and the symmetry of the 3s electron. The dipole moment 
operator transforms according to irreducible representation A'l of the Z^s/j group, since the 
dipole operator is proportional to the spherical harmonic Yim{0' , (p') . As a result, the final 
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state of the electron-ion complex must have A" total symmetry. Since we only consider final 
electronic states with / = 1 (the symmetry is A", when the electron is at large distances 
from the ion), the final symmetry of the ion should be AL. 

As in to the previous theoretical treatments p, lO, III, in order to evaluate the dipole 
transition moments from the initial states in Eqs. ()28j] and ^I^i . we use the 3s molecular 
potential surface of H3 calculated by Nager and Jungen \^ , in a Coulomb approximation. 
Vibrational wave functions of the {00"} and {10°} initial states were calculated using the 
adiabatic hyperspherical approach as described above. 

B. Scattering wave function of e +H^: all asymptotic channels are open 

In the previous section, the scattering matrix for electron-ion collisions have been pre- 
sented. However, for determination of transition dipole moments we need to know not only 
the scattering matrix, but also wave functions of corresponding scattering states. Following 
.e„e.a, defect ..eo.. 013, we wit. a .aUen,. g wave function assuming 

that the electron energy is so large that all possible entrance channels are open. We use the 
same phase conventions for wave functions as in Ref. The wave function having 

out-going wave in the channel i only can be represented as a Ntot component vector, where 
each component "^u/ with i = 1, ■ ■ ■ , Ntot corresponds to the incoming wave in channel i: 

i = l,--- ,Ntot. (30) 

The functions are outgoing/incoming waves in channel i. They are defined in Ref. jl^ . 
The factor ^i{uj) is a part of the total wave function; u includes all degrees of freedom exclud- 
ing the radial one, r. The wave function ()30|) can be considered as a complex conjugation of 
(or time-reversed to) the familiar incoming wave scattering state, having an incoming wave 
only in the channel i'. There are Nfot functions of the type (j30j) and the whole set of Ntot 
wave functions with Ntot components can be considered as a Ntot x Ntot matrix ^. 

C. Scattering wave function of e~+Ht when some channels are closed 



When the energy of the system is low enough such that some asymptotic channels are 
closed to ionization, the total wave function of the system must asymptotically (in the radial 
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coordinate r) vanish in the corresponding channels. Thus, the total wave function differs 
from the one given by Eq. (jHUj) . Let No and Nc = Ntot — No represent the numbers of 
open and closed channels at a given total energy E. In this situation: (i) there are only No 
physically acceptable wave functions of the type (j3(J|) instead of Nfou (ii) these No functions 
are zero at infinity in closed channels. Every acceptable wave function should have the 
following asymptotic behavior [l^ : 



^i;) ^0, i = No + lr-- ,Ntot (31) 



As is well-known, and shown in Refs. 
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12| |. one way to obtain states with this asymptotic 
behavior is to construct linear combinations of states In the matrix form this can be 
written as 

= *B. (32) 

The Ntot X No matrix consists of No vectors, each having Nt„t components. The Ntot x No 



matrix B of the linear transformation is derived in Ref. 
matrix into open and closed subspaces, as 



12| . If we partition the coefficient 




(33) 



the open-channel part Bq is represented by an No x No identity matrix and the closed part 
Be is 

Be = -(SL-e^^^)-'sL (34) 

In the above equation the matrices S^^, and S^.^ are submatrices of 5^ , which is itself 
partitioned as 



(35) 




and (3{E) is a diagonal A''^ x A^,, matrix 
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where Ei refers to a particular ionization threshold i = {vi, f2^}(A^+, K~^). 

After we apply the transformation ()H2|1 . the component \E'^j,'* of the i-th independent wave 
function in the open channel o is given outside the reaction volume by 

o=l,---,iVo, (37) 
where the physical No x scattering matrix 5^^^^* is 

5tP/... = SL-SL(SL-e2^^)-^SL. (38) 
The closed-channel components \E'|.^'' of ^[7^ are determined by 

c= No + l,--- ,Ntot. (39) 

In the above equation, Wc{r) is the Whittaker function and Z is the A^^^ x No matrix of 
closed-channel coefficients (see Eqs. f_2.52-2.54) of Ref. jlj]). In a compact notation, the 
wave function ^^v"'* can be written as 

No . 

^(-) = E'^°(^)-7f {foir)So,' - fo-{r)Slfy') + 

o=l ^ 

Ntot 



^ ^,{uj)W,ir)Z,,,. (40) 



c=No+l 

D. Dipole moments of transitions from the initial bound state of H3 to scattering 
states 

We need to evaluate No dipole transition moments from a fixed initial state, "^ini = 
[Ss^A'^, {00°}(1,0)] or [Ss^A'^, {10°}(1,0)] into all No final states Each such moment 

is represented as 

df = {^-^\e-f\^,m), (41) 

where e*is a unitary vector of laser light polarization. The initial state is closed for autoion- 
ization and, therefore, is represented as \l/i„j = $j„j(ct;)VFj„j(r)z/j„^''^. The factor I'in/'^ is due 
to the unity normalization of the initial bound state: The Whittaker function itself vanishes 
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at infinity, but is chosen to have an energy-normahzed amphtude at small r. Therefore, the 
dipole moment df is written as 



0=1 



Ntot 



j=No+l 



(42) 



The dipole moment df above depends strongly on energy and, therefore, it must be 
calculated at a fine energy mesh. Calculation of all terms in Eq. 021 is computationally 
expensive. However, inspecting the terms in Eq. 021 we notice, that each term can be 
represented as a product of two factors: one factor depends strongly on energy, another 
factor is weakly energy- dependent. Briefly, the two sums above can be combined in one 
single sum of the form: 



Ntot 

df = J2^^^if' 

j=0 

where each term of the sum is represented as product of two factors dj and ■Ujj 



dj 







{^j{u;)W,{r)\e-f\^U^)WUrKnT), if J = iV^ + 1 . . . iV, 



(43) 



(44) 



Hot 



1, if j = 

_SjPhys^ ifj = l...iV„ 

Z,j, iij = No + l...Nt,t 



(45) 



Notice that in Eq. (jl^ the summation starts at j = but not at j = 1 as in Eq. (021) • This 
is because we represent the term with j = / in (Eq. as a sum of two terms in Eq. ()43p. 
with j = (corresponds to f^{r)6jj) and with j = f (corresponds to f~{r)Sj^^^^). This is 
necessary in order to separate energy dependent factors from energy independent ones. 

The partitioning given by Eq. (j43p allows us to significantly reduce the calculation time, 
because over the small range considered in this calculation, it is adequate to evaluate dj only 
once for all energies. However, dj should be calculated at every energy point. 
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E. Transition dipole moment: Integration over all degrees of freedom 



Calculation of Ntot dipole moments of Eq. ()42|) implies integration over all degrees of 
freedom. In practice, the integration is accomplished in several steps. First, we integrate 
over radial coordinate r, calculating geometry dependent elements and nic- 

m^(Q) = ±— =/ ft{r)WUr>iTrdT, (46) 



In order to avoid strong energy-dependence of these elements, we do not include factors 
S]^!^^'^ and Zci' at this stage of the calculation. Therefore, the elements and nic are 
essentially independent of the photoionization energy. 

In Eq. the Whittaker function Wini{r) is calculated for / = and the effective 

quantum number i^^n/^{Q) that depends of configuration and linked to the 3s potential of 
Ha as 

^^ni{Q) = ^ (47) 

where V'^i^Q) is the ionic potential and V^'^{Q) is the potential of the 3s state. Quantum 
defects needed to calculate functions Wdr) and ff{r) of final states are obtained from 
diagonalization of the scattering matrix S. The integral of Eq. is calculated numerically. 

When evaluating Eq. ()46|) . we must choose the principal quantum number n for final 
states (n is always 3 for the initial state). For a given energy E of the photo- ionizing neutral 
molecule, "n./ is determined individually for every final state / using ionization threshold 
energy E~j of the channel /. \{ E approaches closely to the threshold E^ , we take a large 
but finite n, typically n = 20. Such approach is justified by the fact that the overlap in 
Eq. fl4(i|l depends weakly on n/ if n/ is significantly larger than = 3. In principal, this 
procedure can be used for each photoionization energy E. Since the typical theoretical 
spectrum is calculated for typically more than 10^ energy points, this makes the evaluation 
of Eq. at all energies E expensive. To reduce the calculation time, we have adopted 
fixed values of Uf and, therefore, fixed overlaps in Eq. ()4(ij) for several energies E. Since 
rrij (Q) in Eq. ()4fi|l varies slowly with E, this approach provides much more rapid and 
sufficiently accurate values of the matrix elements in Eq. 

The next step in our evaluation of the dipole moments of Eq. ()42|) is the integra- 
tion over the vibrational coordinates Q. The total vibrational function in the adiabatic 
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approximation is represented as a product of hyper-radial and hyperangular components 
= ipv^{R)(^v2i2i(^A)- Knowing the wave functions '^"'■^{Q) and the Q-dependent 
matrix moments of Eq. (|46|) . we can calculate the desired dipole matrix elements between 
initial and final vibrational states as follows: 

R {{vl,v2'^},- {vl,v2%^,) = JdQ ^'^Xv2^2y{Q) (Q) vl/«^(Q). (48) 

where mj is mf or rric. In the same manner as for the scattering matrix, we evaluate this 
integral in two steps; first in the space of hyperangles, then in the space of the hyper-radius. 

The next step is to evaluate the integration over angular coordinates in ()44j) . To do this, 
we write explicitly all quantum numbers of the final and initial wave functions, which is 
given by Eq. ()24|) . and we write the dipole moment dj as 

= (iV+, K+; Nj, m,{vu v'i}j\e- f| Ar+., K+., AT^^,, mini{vuv'^'}ini)- (49) 

This expression can be represented in a form that is more suitable for our calculations. After 
some angular momentum algebra (for more details see the appendix) we derive the following 
formula for the dipole moment dj. 



^ ~ ^j2Nj + 1 

x$:(-i)^-^^<(t^.A.<A^to^({^i'^^n.; (50) 

A. 

It gives the dipoles moments dj in terms of vibrational matrix elements of Eq. EHl which 
are calculated numerically. 

In applying the present treatment to the interpretation of the photoionization experi- 
ments by Bordas et al. [4t] and by Mistrik et al. [5], we consider only two different initial 
states specified by the set of quantum numbers i = [{vi, V2}inikni-^ini]- For the final state, 
[{vi,vl^}jljAj], only the index Ij = 1 is fixed — we consider only p-wave final states. 



F. Final theoretical photoionization spectrum 

After the matrix elements dj have been obtained, the dipole transition moments df are 
calculated using Eq. As we mentioned before, dj can be calculated once and used for all 
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photoionization energies, provided the energy range of the entire calculated photoabsorption 
spectrum is not too extensive. However, the coefficients vjf must be recalculated at every 
final state energy of the theoretical photoionization spectrum. Having calculated the dipole 
transition moments dj, the total oscillator strength into the open ionization channels is then 
given by 



^-s-Ew'. (51) 

/=i 



where oo is the frequency of the laser light. 

G. Quantum defect parameters used in the calculation 

In order to construct the scattering matrix and dipole moment vector, we use the quantum 
defect parameters 6, \, fi\=o{Q), and fiA=±i determined by Mistrik et al. 5] from accurate 
ah initio calculations of potential-energy surfaces of H3. As mentioned in our previous work 
0, yi , the parameters 5, A determined in Ref . should be multiplied by factor — tt in order 
;o correct the convention inconsistency in the definition of the reaction matrix K in Ref. 
5|. The parameters 5, A, yUA=o(Q)5 /^a=±i are slightly different for different Rydberg states 
5|. In our calculation we use values obtained for n = 4 Rydberg states: 5 = —n ■ 1090 
cm~^, A = — TT • 12360 cm~^, fi\=±i = 0.395. The quantum defect /iA=o(Q) depends weakly 
on nuclear configuration Q. We use /iA=o(Q) from Ref. j^. The ionization energies of 40 
lowest states are taken from the same reference. 

V. RESULTS AND COMPARISON WITH THE EXPERIMENTS 



In the experiment by Bordas et al. [4!], the initial state of H3 is the state ii = 
[H3 A'2{tot.){00°}{N = 1,K = 0), 3s], the initial state in the experiment by Mistrik et 
al. js] differs only by the symmetric stretch vibrational quantum number Vi, which is the 
singly excited: = [H3 A'2{tot.){10^}{N = 1,K = 0), 3s(e/.)]. The energy difference be- 
tween these two states is 3212.6 cm~^. In the first experiment the energy region around 
the ground rovibrational level of the ion is probed by a tunable laser. The energy difference 
between the state ii and the ground rovibrational statejHj]' A2{ion) {00^} {N~^ = 1, = 0)] 
of the ion is 12867.6 cm~^. In the second experiment [s], the energy region around the state 
with singly excited f 1 = 1 mode [H3 (ion) { 10°} (A^+ = 1,K^ = 0)] of the ion is probed. 
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The energy difference between these two ionic levels, [Hj]' A2{ion) {00^} {N~^ = 1,K^ = 0)] 
and [Rj A'^lion) {10°} {N^^ = 1,/^+ = 0)], is 3176.06 cm'^ In the present treatment the 
energy origin is set to the ground ionic state. 

Figures Q and shows experimental and calculated spectra for the two experiments. 
The overall agreement is quite good. Below we give a more detailed discussion of the 
comparison between the experiments and our calculation. 

In constructing the theoretical spectrum, we have combined the spectra for = and 
N = 2 according to the experimental conditions. Specifically, to compare our theoretical 
results with the experiment by Bordas et ai, we have summed up the separate theoretical 
spectra for = and N = 2. To compare with the experiment by Mistn'k et al. we have 
accounted for a fixed angle of 60° between the linear polarization vectors of two the lasers 
used in the experiment. We used the prescription of Ref. 



5^3, according to which the final 



spectrum ^ is constructed as 



dE 4dE 16 dE' ^ ' 

where ^ and ^ are the spectra calculated for = and N = 2, respectively. 

The energy regions accessible in the experiments have three qualitatively different 
regimes, namely the discrete, Beutler-Fano, and continuum regimes. 

The Beutler-Fano regime arises at energies between two different rotational levels asso- 
ciated with the same vibrational state. For both experiments, this region occurs at final 
state energies between the [{viO°}{N+ = 1,K+ = 0)] and [{viO°}{N+ = 3, K+ = 0)] states 



of the ion, where Vi = for the experiment of Ref. 



^ and fi = 1 for the experiment 



of Ref. lO]. Autoionization in this region usually occurs quite rapidly compared to au- 
toionization in other energy ranges. Consider, for example, an electron is excited into a 
Rydberg state attached to the [{f iO°}(A^+ = 3, = 0)] ionic level, at a total energy above 
{fiO'^}(A^'^ = 1, = 0)]. Since both rotational levels have the same vibrational excitation, 
the corresponding Franck-Condon overlap between the Rydberg state and a continuum state 
of the {viO°}{N^ = 1,K^ = 0)] level is very favorable, which generates a large autoioniza- 
tion width, as is evident in the Beutler-Fano regions of both experiments. Figure El shows 
a detailed comparison between theory and the experiment ^ for the Beutler-Fano region, 
and Fig. IHl presents a detailed comparison with the second experiment j^. The agreement 
between theory and experiment is good, which is evidence supporting the approximations 
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we have adopted in our theoretical description. 

The continuum regime is situated above the corresponding Beutler-Fano energy range. 
Generally, autoionization is much slower in such regions. Figure [7| compares the calculated 
spectrum with the experiment of Ref . ^ . We draw attention to two broad resonances around 
740 cm^^ and 950 cm^^. Not only do these nicely reproduce the experimental spectrum, but, 
importantly, they are caused by Jahn- Teller coupling. However, in the previous experimental 
and theoretical studies, these features were ignored, probably because they were construed 
as noise. The present calculation suggests that these are real resonances, broadened by a 
strong interaction between rotational and vibrational degrees of freedom. Figure |H1 compares 
the present calculation with the results of Ref. In this figure we would like to note two 
other features. The experimental data is rather noisy, but an inspection of the calculated 
and experimental spectra suggests that the calculated resonances appearing around 3770 
cm~^ and around 3850 cm~^ correspond to broad observed resonances. Again, the large 
widths of the two resonances are caused by Jahn- Teller coupling. 

The discrete regime is the energy range where autoionization of Rydberg states is ener- 
getically forbidden. However, we will follow the convention proposed in Ref. and will call 
the region below [H^ A'2{ion) {00^} {N^ = 1, = 0)] as discrete too when we refer to ^2 as 
the initial state. This convention is justified by the fact that autoionization of such states 
is slow, owing to an unfavorable Franck-Condon overlap of these states with the ground 
ionic state. Figure El shows a comparison of experimental and calculated spectra for the 
discrete spectrum of the i2 initial state. 

Although the overall agreement between theory and experiment is still good overall, 
there are some resonances in the experimental spectrum that absent from the calculated 
spectrum. Some of these missing resonances could conceivably be caused by an influence of 
the d electronic wave. This possibility was demonstrated by Mistrik js^. For example, an 
additional resonance in the experimental spectrum around 2970 cm~^ appears to be caused 
by an nda Rydberg electron. 

VI. CONCLUSION 

In the present study, we propose an updated theoretical method for treatment of photoion- 
ization in the H3 molecule. The main engine of the method, MQDT including the Jahn- Teller 
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coupling and the rovibrational frame transformation, is the same as in studies by Stephens 



and Greene 



6, 



in the treatment of the 



7[. However, we have improved beyond Refs. 
symmetry issues associated with different degrees of freedom. We have proposed a new and 
more efficient symmetrization procedure that accounts for the symmetries of nuclear spin as 
well as the rotational and vibrational parts of the total wave function. Another improvement 
is that we include into the description of H3 photoionization the possibility that the system 
may break apart by dissociation in addition to ionization. To our knowledge, this is the 
first time that photodissociation has been included in its competition with photoionization 
for richly resonant Rydberg spectrum of a triatomic molecule. Although in this system, 
predisso elation apparently does not play an important role for the overall spectra under 
consideration, it must be important, at least, for some Rydberg states of H3. In fact, the 
importance of the dissociation channel was recently demonstrated experimentally . In a 
future study, hopefully, we will investigate theoretically the predisso elation of such Rydberg 
states. The inclusion of predissociation could also be important in studies of other triatomic 
molecules. Note that as in Ref. 0, |^, we have corrected the earlier inconsistency in the 



reaction matrix convention employed in Refs. M, l2l, 12^ which results in a multiplication of 
the Jahn- Teller parameters A and 6 employed in those references by the factor of — tt. In the 
present study, we have also improved the dipole moment calculation, and reformulated it in 
a form that should be particularly accurate for the photoionization of a Rydberg molecular 
initial state. In contrast to most previous studies of H3 photoionization, 0, Q] the present 
treatment makes use of the adiabatic hyperspherical approach for the vibrational motion 
of the nuclei in the Hjj' ion and in the H3 initial state being photoionized. Owing to the 
adiabatic approximation, this approach should be a priori less accurate than the one used 
in Refs. Q], but the loss of accuracy due to this approximation seems to be small. 

The possibility of constructing a unified theoretical description of H3 photoionization and 
Ht dissociative recombination is an attractive feature of this approach. 

nn 

We have obtained good agreement with both photoionization experiments |^ y] and with 
many of the spectra obtained in previous theoretical studies j^, 0, Q]- In some regions the 
agreement with the experiments is even better in the present study than it was in Refs. 
0, y, Q]- We attribute this to the correction of the aforementioned error in determination 
of the Jahn- Teller coupling parameters A and 6. Although we have been able to reproduce 
most of the observed resonance features, both in position and in shape, there remain several 
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features in the experimental spectra that are not described by our treatment. One possible 
explanation of the discrepancy between theory and experiment is the influence of nd Rydberg 



states, as proposed by Mistn'k [39[, because we have not taken these states into consideration. 

In conclusion, we would reiterate that we have developed a new theoretical method that 
can be applied to a unified treatment of photoionization and dissociative recombination for 
molecules of the D^^h symmetry group. Our application to the H3 molecule shows good gen- 
eral agreement with existing experimental observations of H3 photoionization. The method 
has also been shown in our previous study 0, to give good agreement with measurements 
of Hg" dissociative recombination. 

VII. APPENDIX 

In the appendix we provide a detailed derivation of Eq. EHl starting from Eq. 0^ 
Suppose the laser light is linearly polarized in the Z direction in LS. The covariant 

spherical tensor components of the polarization vector are then (e+i, eo, e_i) = (0, 1, 0). The 

vector r in spherical coordinates is js^ 

47r 

y|r|Fi,^(^,0) with/i = 0,±l. (53) 
Then, the scalar product e ■ r in LS is 

(e-rO,.5 = yf |r|yi,o(^,</>). (54) 
The expression of this product in MS is obtained by an appropriate coordinate rotation. 




(e ■ r)MS 



^\r\^Y,^^,{e\<P') [Pi^,(a,A7)]*. (55) 

Accounting Eqs. ()25|) and (j2Sl); the expression for the amplitude of the dipole transition of 
Eq. becomes 

A' 

^ |-^jni5 ^tni'i ^inii ^imf'^lj V2}ini) = 

V ^ A',A, 

X [Pi^,(a/57)]* 7^(iV,„„ m,„„ K,^ aPl)Yu„.K^XO' . ^)) ^,,e','CB:^,N.„.,K.„. x 

xi-iy-^-''-^R{{vi,v'i}j; {vi,vl,'}ini) . (56) 
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If the initial electronic state is 3s, then Yi.^.i^-^- = Yq^q = l/v47r. The integral over electronic 
angles 9' , yj' is then trivial, and the relevant angular matrix element is 



6, 



(57) 



The integral over a,/3, and 7 angles is evaluated using the addition theorem for Wigner 
functions and their normalization properties 3^ : 

\^mj ,Kj I -^0,A' I ^mi„i,Ki„i / 



\2Nj + l){2Nj + l) 



Stt 



2^2 



a/37 



'{2Nj + l){2Nini + l] 



(8n 



2\2 



\^mj,Kj\-^mini,Kini+K'^lfi;Ni 



N' 



a/37 



'2Nini + 1 



^'-^l,0;Af™,m,„i<-'l,A';Ar,„i,K,„,('Afi,Af'f'm,,m,„i('it,,K,„,+A' - 10»J 



-,NT,Ki„i+K' 



2N,+ 



The matrix element then becomes 



do 



1 2M., + 1 



X 



^lj-Aj;Nj,Kj'~^l,0-,Ni„i,mi„i'~^l,Aj;Ni„i,Ki„i"rnj,mi„iOKj,Ki„i+Aj^ X 

A. 

(-1)1-^^7? 4^},; {vuv'ihr,,). (59) 

Using the fact that Ni^i = N^^ = 1, Ki^i = Kf^^, rriini = mf^^ = mj, Kj = Ki^i + Aj 
and also the requirement K^- = 0, which arises from symmetry restrictions on the initial 
rovibrational state, the last relation simplifies in this case to: 



1 ^Nj ,mini 

'-'l,0;l,mi„i ^ 



^/2Nj + 1 



X 



(60) 



The last formula is Eq. EHl 



28 



Acknowledgments 



This work has been supported in part by NSF, by the DOE Office of Science, and by 
an allocation of NERSC supercomputing resources. The authors are grateful to H. Helm, 
M. Larsson, J. Stephens, and I. Mistrik for fruitful discussions and to E. Hamilton and B. 
Esry for assistance. 



[1] H. Helm, Phys. Rev. Lett. 56, 42 (1986). 
[2] H. Helm, Phys. Rev. A 38, 3425 (1988). 

[3] A. Dohdy, W. Ketterly H. P. Messmer, and H. Walthcr, Chem. Phys. Lett. 151, 133 (1988). 
[4] M. C. Bordas, L. J. Lembo, and H. Helm, Phys. Rev. A 44, 1817 (1991). 
[5] L Mistrik, R. Reichle, U. Miiller, H. Helm, M. Jungen, and J. S. Stephens, Phys. Rev. A 61, 
033410 (2000). 

[6] J. A. Stephens and C. H. Greene, Phys. Rev. Lett. 72, 1624 (1994). 
[7] J. A. Stephens and C. H. Greene, J. Chem. Phys. 102, 1579 (1995). 
[8] V. Kokoouline and G. H. Greene, Phys. Rev. Lett. 90, 133201 (2003). 
[9] V. Kokoouline and G. H. Greene, Phys. Rev. A 68, 012703 (2003). 
[10] M. J. Seaton, Rep. Prog. Phys. 46, 167 (1983). 

[11] U. Fano and A. R. P. Rau, Atomic Collisions and Spectra (Academic Press, Orlando, Florida, 
1986). 

[12] M. Aymar, G. H. Greene, and E. Luc-Koenig, Rev. Mod. Phys. 68, 1015 (1996). 
[13] Gh. Jungen Molecular Applications of Quantum Defect Theory, (Institute of Physics Publish- 
ing, Bristol, U.K., 1996). 
[14] O. L Tolstikhin, V. N. Ostrovsky and H. Nakamura, Phys. Rev. Lett. 79, 2026 (1997). 
[15] O. I. Tolstikhin, V. N. Ostrovsky and H. Nakamura, Phys. Rev. A 58, 2077 (1998). 
[16] E. L. Hamilton and G. H. Greene, Phys. Rev. Lett. 89, 263003 (2002). 

[17] B. J. McGah, A. J. Huneycutt, R. J. Saykally T. R. Geballe, N. Djuric, G. H. Dunn, J. Sema- 
niak, O. Novotny, A. Al-Khahli, A. Ehlerding, F. HeUberg, S. Kalhori, A. Neau, R. Thomas, 
F. Osterdahl, and M. Larsson, Nature (London) 422, 500 (2003). 

[18] T. Tanabe, K. Chida, T. Watanabe, Y. Arakaki, H. Takagi, I. Katayama, Y. Haruyama, 



29 



M. Saito, I. Nomura, T. Honma, K. Noda, K. Hoson, in Dissociative Recombination: Theory, 

Experiment and Applications IV, edited by M. Larsson, J. B. A. Mitchell, I. F. Schneider, 

(World Scientific, Singapore, 2000) pl70. 
[19] M. J. Jensen, H. B. Pcdersen, C. P. Safvan, K. Seiersen, X. Urbain, and L. H. Andersen, Phys. 

Rev. A 63, 052701 (2001). 
[20] Y. Zhou, C. D. Lin, and J. Shertzer, J. Phys. B 26, 3937 (1993). 
[21] C. D. Lin, Phys. Rep. 257, 1 (1995). 

[22] B. D. Esry, C. D. Lin, and C. H. Greene, Phys. Rev. A 54, 394 (1996). 
[23] V. Kokooulinc, C. H. Greene, and B. D. Esry, Nature 412, 891 (2001). 

[24] H. C. Longuet-Higgins, in Advances in Spectroscopy, (Interscience, New York, 1961). vol. II 
p. 429. 

[25] A. Staib, W. Domcke, and A. L. Sobolewski Z. Phys. D 16, 49 (1990). 
[26] A. Staib and W. Domcke, Z. Phys. D 16, 275 (1990). 

[27] V. Spirko, P. Jensen, P. R. Bunker, A. Cejchan, J. Mol. Spectrosc. 112, 183 (1985). 
[28] P. Bunker and P. Jensen Molecular Symmetry and Spectroscopy, (NRG Research Press, Ot- 
tawa, Ganada, 1998). 
[29] J. K. G. Watson, Phil. Trans. R. Soc. Lond. A 358, 2371 (2000). 

[30] R. Rcichlc, Ph.D. Thesis, Freiburg (2002). 

[31] W. Gencek, J. Rychlewski, R. Jaquet, and W. Kutzelnigg, J. Ghem. Phys. 108, 2831 (1998). 

[32] R. Jaquet, W. Gencek, W. Kutzelnigg, and J. Rychlewski, J. Ghem. Phys. 108, 2837 (1998). 

[33] G. H. Greene and Gh. Jungen, Adv. At. Mol. Phys. 21, 51 (1985). 

[34] E. S. Ghang and U. Fano, Phys Rev. A 6, 173 (1972). 

[35] G. M. Lindsay, B. J. McCall, J. Mol. Spectrosc. 210, 60 (2001). 

[36] B. J. McGall, Ph.D. Thesis, University of Ghicago (2001). 

[37] G. Nager and M. Jungen, Ghem. Phys. 70, 189 (1982). 

[38] D. A. Varshalovich, A. N. Moskalev, and V. K. Khersonskii, Quantum Theory of Angular 
Momentum, (World Scientific, Singapour, 1988). 

[39] I. Mistrfk, Ph. D Thesis, Bratislava (2001). 

[40] L Mistrfk, R. Reichle, H. Helm, and U. Miiller, Phys. Rev. A 63, 042711 (2001). 



30 



(23) 



(13) 



- . 



.4 .4 

. . .-^ ^ 7 4 r ^» 3. 




1 2 



^ I « 

(12); 




(12)T 



FIG. 1: The figure demonstrates how the two-dimensional space of the hyperangles 6 and (p 
reproduces all possible shapes of a triangle consisting of three nuclei. The size of the system is 
controlled by the hyper-radius, which is not shown here. Hyperangle ip starts at point 99 = and 
goes counter clockwise from to 2it. The hyperangle ^ is at the symmetric configuration and it 
changes to 6 = ti/2 where it represents linear configurations. The figure also shows symmetry axes: 
Every operation in the group C^v involves only the hyperangle ip. For example, the operation (12) 
transforms the triangle Ta into the triangle T;,, which corresponds to a reflection around the line 
(p = 7r/2...37r/2. 
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FIG. 2: Lowest adiabatic hyperspherical potentials of the ion plotted as functions of hyper- 
radius. Potentials of the Ai symmetry are plotted with solid lines, A2 curves are represented 
with dot-dashed lines, and E curves are represented with dashed lines. Every E curve is doubly 
degenerate. For several lowest curves we also specify approximate quantum numbers V2 and l2- For 
excited hyperspherical states such quantum numbers can not be defined since the three-dimensional 
ionic potential is strongly anharmonic at large energies. As is clear in the figure the potential fi4^ 
is already strongly shifted from its partners viA"^ and vi4^ toward to the next family of states uiS'^. 
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FIG. 3: Comparison of calculated photoionization spectra (upper and middle panel) and the 
experimental (lower panel) spectrum obtained by Bordas et al. J]. In this experiment, the initial 
state of H3 is [Ss^Ag, {00°}(1, 0)]. The upper panel shows theoretical spectra calculated separately 
for = (black line) and N = 2 (gray line). The middle panel shows the sum result of the two 
spectra convolved with the experimental resolution width 0.15 cm~^. Therefore, the spectrum in 
the middle panel should be compared with the experimental spectrum. Here and in all following 
figures, the energy is related to the ground rovibrational level of H|^(742), {00'^}(1,0). 
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FIG. 4: Comparison of the theoretical (upper panel) spectrum with the observed one (lower 
panel) from experiment of Ref. ^. In this experiment, the initial state differs from the experiment 
by Bordas et ai: It is [35^^42, {lO'^Kl, 0)]. In constructing the theoretical spectrum, we have 
combined the spectra for = and N = 2 according to experimental conditions (see text). The 
final spectrum is convolved with the experimental resolution width of 0.2 cm^^. 
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FIG. 5: Comparison between theoretical (upper panel) and experimental (lower panel) results: 
the Beutler-Fano part of the spectrum shown in Fig. |31 This energy region corresponds to energies 
between {00'^}A^+ = 3 and {OO^jA^^ = 1 ionic rotational levels. The series of wide resonances 
are due to the fast rotational autoionization of states with = 3 to an open continuum of the 
= 1 rotational level. Interlopers at 70 cm~^ and 150 cm~^ are two examples of resonances 
playing an important role in DR of H^. 
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FIG. 6: Theoretical (upper panel) and experimental (lower panel, from Ref. 
spectra. Theory reproduces the most of experimental features. 



Beutler-Fano 
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FIG. 7: A part of the theoretical (dark hne) and the experimental (gray line) spectra (see Fig. I^J 
corresponding to the H3 continuum above the {00'^}(30) rovibrational level of Hj^. Relatively wide 
resonances at 740 cm~^ and 950 cm~^ are produced by the Jahn- Teller coupling between electronic 
and vibrational motion of H3. Although the overall agreement between theory and experiment is 



better than in previous theoretical studies 
are not reproduced in the present approach 



several experimental features still exist that 
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FIG. 8: Theoretical (dark line) and experimental (gray line) spectra (see Fig. 0} corresponding to 
the H3 continuum above the {10''}(30) rovibrational level of H^. Similar to Fig. [71 the agreement 
between theory and experiment is better than in the previous theoretical study [sj, but there are 
several uninterpreted experimental resonances. 
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FIG. 9: Theoretical (upper panel) and experimental (lower panel) |5( spectra of H3 in the quasi- 
discrete energy region below the {10'^}(10) rovibrational level of H^. The initial state for the 
dipole transition (H3 photoexcitation) is [Ss^^dg, {10''}(1, 0)]. This region is referred as quasi- 
discrete since states, populated starting from [Ss^Ag, {10'^}(1, 0)], are only weakly interacting with 
lower ionic states with v = {00°}. 
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FIG. 10: The figure demonstrates a correlation between photoionization (upper panel) and dis- 



sociative recombination (lower panel) 



spectra obtained theoretically and in the experiment. The 



n 



upper panel shows both experimental ^ and theoretical spectra, similar to Fig. |31 The lower panel 
gives an experimental DR rate obtained in Ref. which is shown as circles. Grey line in lower 
panel gives theoretical probability of DR process, calculated for the total A'2 symmetry of H^. 
This symmetry corresponds to the total symmetry of H3 in the both photoionization experiments 
considered in this study. 
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